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ABSTRACT 

We report the results of high spectral and angular resolution infrared obser- 
vations of the multiple system GL 569 A and B that were intended to measure 
the dynamical masses of the brown dwarf binary believed to comprise GL 569 B. 
Our analysis did not yield this result but, instead, revealed two surprises. First, 
at age ~ 100 Myr, the system is younger than had been reported earlier. Second, 
our spectroscopic and photometric results provide support for earlier indications 
that GL 569 B is actually a hierarchical brown dwarf triple rather than a binary. 
Our results suggest that the three components of GL 569 B have roughly equal 
mass, ~ 0.04 M© 

Subject headings: stars: binaries: visual — binaries: spectroscopic — stars: low- 
mass, brown dwarfs — stars: pre-main sequence — techniques: high angular 
resolution — techniques: radial velocities 



1. Introduction 

Accurate and precise measurements of the masses of brown dwarfs are necessary to 
validate the theoretical calculations of their evolution and to determine their location in 
the HR diagram (Burrows et al. 1997; Baraffe et al. 2003). Astrometric and spectroscopic 
observations of binaries provide the means to accomplish this dynamically. When the dis- 
tance to the binary is known, the astrometric orbit of an angularly resolved binary (hence a 
"visual binary", VB) determines the total mass, Mtot- A double-lined spectroscopic binary 
(SB2) provides the mass ratio. When a binary can be observed as both a VB and SB2, the 
component masses can be determined. Since the VB also yields the orbital inclination, the 
VB and SB2 measurements together provide an independent determination of the binary's 
distance. In the special case of an eclipsing binary, the orbital inclination is ~ 7r/2 and the 
SB2 observations alone suffice to determine the component masses. 
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Although brown dwarf binaries are not rare (Burgasser et al. 2005), precise dynamical 
masses are known only for the brown dwarf components of the eclipsing binary in Orion 
discovered by Stassun et al. (2005). Their masses, radii, and effective temperatures agree 
reasonably well with the models, but surprisingly, the effective temperature of the primary 
is found to be cooler than that of the secondary. Zapatcro Osorio et al. (2004), Bouy et al. 
(2004), and Pravdo et al. (2005) have reported brown dwarf masses based on an astrometri- 
cally determined M^ot and photometry, but reliable masses determined dynamically are not 
yet available. Zapatero Osorio et al. (2004) obtained angularly resolved spectra of GL 569 Ba 
and Bb but their measurements were insufficient to yield dynamical masses. Because only 
two brown dwarf masses are known, and they raise a puzzle that requires confirmation, it 
is desirable to enlarge the sample of dynamically determined masses and to cover as large a 
range of mass and age as possible. 

GL 569 is a multiple system comprised of GL 569 A, a M2.5V star (Henry & Kirkpatrick 
1990) of mass M = 0.35 ± 0.03 Mq (Gorlova et al. 2003) at the HIPPARCOS distance 
9.81 ± 0.16 pc, and GL 569 B, thought to be a binary brown dwarf with period ~ 876 days 
(Zapatero Osorio et al. 2004). We began a study of GL 569 B in 2002 because the astrometric 
orbit is well determined and the components, Ba and Bb, are bright enough to be observed 
at high spectral resolution in the near-IR. Thus, GL 569 B is well suited for the fusion of VB 
and SB2 techniques to measure the mass of its components. This paper presents our results. 
We describe our observations and data reduction in §2, the results for GL 569 A and B in §3 
and §4, respectively, and discuss our results in §5. We present a summary and suggestions 
for future work in §6. 



2. Observations and Data Reduction 

All observations were obtained at the W. M. Keck Observatory using the facility near- 
infrared spectrometer NIRSPEC (McLean et al. 2000) and near-IR camera, NIRC-2. We 
used NIRSPEC in its high spectral resolution mode centered at 1.55 //m, with and without 
the adaptive optics (AO) system (Wizinowich et al. 2000). The NIRC-2 observations in the 
H- and X-bands were made using AO. Table 1 provides a log of the observations. 

We used a 2-pixel wide slit for all NIRSPEC observations which provided spectral res- 
olution ~34,000 with AO and ~31,000 without it. With AO, the diffraction resolution was 
~ 40 mas. This enabled us to measure the angularly resolved spectra of GL 569 Ba and 
Bb simultaneously with both objects on the appropriately rotated spectrometer slit. All 
observations were taken as a series of 240 — 300 s integrations with the targets shifted in 
an ABBA pattern between two positions separated by ~10" on the slit without AO, and 
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~1" with AO. Sets of 10 flats and 10 darks were median filtered and the results, along with 
individual AB pairs of target data, were input into the REDSPEC code for rectification, 
flat-fielding, and spectral extraction^. The resultant spectra for a single object were then 
averaged. We worked with NIRSPEC orders 48 and 49, centered at ~1.59 and ~1.56;um, 
respectively, because these are the most free of telluric absorption of the 9 orders available 
at the grating setting we used (see Bender et al. 2005, Figure 1). Spectral lamp line inte- 
grations were taken to obtain the dispersion solution for the AO spectra. For the non-AO 
observations, the night sky OH lines were used for this purpose. 

We measured radial velocities of the targets with the two-dimensional correlation algo- 
rithm TODCOR (Zucker & Mazeh 1994), following procedures we have described previously 
(e.g., Mazeh et al. 2002). TODCOR requires template spectra for the primary and secondary 
components. When we observed and reduced NIRSPEC template spectra for the first time 
(Prato et al. 2002), we compared their radial velocities to similar measurements of the same 
stars by Duquennoy & Mayor (1991) and Latham (1985) in order to determine the uncer- 
tainty, ±1.0 km s~^, in our velocities. More recently, Nidever et al. (2002) have published 
very precise radial velocities of 844 FGKM stars, including several whose IR spectra we use 
as templates. We used the Nidever et al. (2002) measurements to correct the velocities of 
stars earlier than M4 in our suite of templates and used the corrected values to bootstrap 
improvements to templates of later spectral type. We hst the revised values in Table 2; for 
the stars in common with Prato et al. (2002) all the revised velocities are within 2.2 km s~^ 
of the values given there, and most agree to a few tenths of a km s~^. 

The NIRC-2 images have a scale of 0'.'00942 per pixel and a field of view of ~9'.'6, sufficient 
to include both GL 569 A and B simultaneously. We took data in the H- and X-bands in a 
five-point dither pattern. Integration times for each frame varied from 0.18 — 1.00 s. After 
fiat- fielding the individual images with the dark-subtracted fiat, we expanded the the pixel 
scale by a factor of 10 to improve the measurement precision of the relative positions and 
photometry for GL 569 A, Ba, and Bb. For the February, 2005 data, we used GL 569 A to 
ccntroid the images for shifting and adding multiple exposures and to create an instrumental 
point-spread function (PSF) in the H- and fC-bands. GL 569 A was saturated in the images 
obtained on December 2004 so we could use these data only for astrometry and relative 
photometry of GL 569 Ba and Bb. 



See http://www2.keck.hawaii.edu/inst/nirspec/redspec/index.html 
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3. Results: GL 569 A 

The proper motion of GL 569 A is /i« = 275.95 and fis = -122.12 mas yr^^ {HIPPAR- 
COS). Thus, in the ~20 years since the discovery by Forrest et al. (1988) of GL 569 B, the A 
component has moved ~6" on the sky. Our February, 2005 NIRC-2 H- and i^-band images 
showed the photocenter of GL 569 B 4'/18±0'/21 North, and 2'.'19±0'/14 East with respect to 
GL 569 A. This offset of B with respect to A hes ~1" from that measured by Forrest et aL 
(1988) (Table 3, Figure 1) and confirms their finding that A and B form a common proper 
motion system. Our measurement, taken together with the positions measured by Martin 
et al. (2000) and Lane et al. (2001), clearly indicates the orbital motion of GL 569 B with 
respect to A (Figure 1). 

Using spectroscopy in visible light, Henry & Kirkpatrick (1990) showed that GL 569 A's 
spectral type is M2.5 V. To measure the radial velocity of GL 569 A and, in the process, to 
estimate its spectral type, we analyzed the order 49 spectra obtained on UT 20 April 2003, 
UT 14 June, 2003, and UT 25 May 2004. The spectrum of the template star GL 436, spectral 
type M2.5 V (Prato et al. 2002), produced the highest correlation with GL 569 A's spectrum, 
in agreement with Henry & Kirkpatrick's result. The difference between the average of our 
radial velocity measurements, —7.7 ± 0.6 km s~^ (Table 4), and the center-of-mass velocity 
of the Ba— Bb pair (§4.3) is 0.8 ± 0.9 km s~^. The velocity difference arising in the orbital 
motion of B with respect to A is expected to be small (§5) and to lie below the significance 
limit of our measurements. 

The apparent magnitudes of GL 569 A at F and H are 10.200 ± 0.004 mag and 5.99 ± 
0.02 mag, respectively {HIPPARCOS, 2MASS), corresponding to absolute magnitudes My = 
10.24 mag and Mh — 6.03 mag at distance 9.81 pc. Using near-IR surface gravity indicators 
and model stellar spectra, Gorlova et al. (2003) estimated that GL 569 A's mass is 0.35 ± 
0.03 M0. Figure 2 shows My and Mh vs mass diagrams using these values. Surprisingly, 
these mass-luminosity diagrams indicate an age of 20—40 Myr for GL 569 A, much younger 
than the several hundred Myr estimates of the ages of the brown dwarfs Ba and Bb (Lane et 
al. 2001; Zapatero Osorio et al. 2004). It is possible, however, that the theoretical isochrones 
underestimate the age of a star in this range of mass and age (1. Baraffe, priv. comm.). 
To obtain an independent estimate of the age, we followed K. Luhman's suggestion (priv. 
comm. 2005) and placed GL 569 A on a (V— K) color vs. absolute K magnitude diagram 
(CMD) for the Pleiades cluster, thought to be 100-125 Myr old (Meynet et al. 1993; Stauffer 
et al. 1998). Figure 3 shows that GL 569 A lies on the Pleiades main-sequence, indicating 
that it is about the same age. GL 569 A could be older than it appears in Figures 2 or 3 if it 
were fainter and its apparent magnitude represented that of two stars in a binary. However, 
its spectrum in both the visible and IR is that of an M 2.5 with no evidence of either a 



- 5 - 



companion, or variability of its velocity on the three occasions that we measured it in the 
IR. Its H and K-band images also show no evidence of a companion. We conclude that 1) 
GL 569 A is a single star, 2) its age is 100-125 Myr, and 3) GL 569 A seems to be younger 
than the age previously estimated for its brown dwarf companions. 

4. Results: GL 569 B 

4.1. Photometry, Astrometry, and Orbital Parameters of GL 569 B as a 

Visual Binary 

Figure 4 shows the H- and i^-band composite images of GL 569 Ba-Bb from the UT 
25 February 2005 NIRC-2 observations. We modeled the brown dwarf binary as the sum of 
two scaled PSFs; the single star PSF was provided by the simultaneously obtained image of 
GL 569 A, which had FWHM = 41 mas at H and 53 mas at K. The AO correction achieved 
a Strehl ratio of ~0.3 at H and ~0.4 at K. We determined the positions and amplitudes of 
Ba-Bb relative to each other, and to GL 569 A, by least squares fitting. We used the 2MASS 
magnitudes of GL 569 A at if and 5.99 ± 0.02 mag and 5.77 ± 0.02 mag, respectively, 
to determine the magnitudes of Ba and Bb. The astrometric and photometric results are 
given in Table 5. Because we could not use GL 569 A as a PSF (§3) in the December, 2004 
data, we formed a composite Ba-Bb PSF by cross-correlating the images and shifting them 
for greatest correlation. We measured the position of Bb relative to Ba by modehng the 
PSF as an eUiptical Gaussian and fitting to the composite image. The astrometric results 
and Bb/Ba flux ratios from the December, 2004 observations are included in Table 5. Ba is 
0.53 ± 0.08 magnitudes brighter than Bb at K, consistent with the measurement of Lane et 
al. (2001), and 0.61 ± 0.07 magnitudes brighter at H. 

Observations of a VB yield P, a("), e, i, Tq, o;, and Q - the orbital period, semi-major axis 
in arcseconds, eccentricity, inclination, time of periastron passage, longitude of periastron, 
and longitude of the ascending node, respectively^. When the distance to the binary is 
known, as it is here, the semi-major axis, in physical units, and period determine the total 
mass of the binary via Kepler's Third Law. Figure 5 shows the apparent orbit derived from 
the astrometry reported by Lane et al. (2001), Zapatero Osorio et al. (2004), and Table 
5. We calculated the orbital parameters using procedures written by Schaefer (2005) and 
described in Schaefer et al. (2003). The new values (Table 6) are in excellent agreement 



and Q. determined for a VB are degenerate by 180°; the ambiguity is removed by observation of the 
system as a spectroscopic binary. 
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with those published by Zapatero Osorio et al. (2004). Our observations extend the time 
over which the orbital motion of the GL 569 B binary has been monitored from 0.85 to 2.3 
orbital periods; as expected, their main effect is to improve the precision of P and Tq. The 
total mass of Ba— Bb that follows from these values is Mtot = 0.125 ± 0.007 Mq. This value 
is consistent with that determined by Zapatero Osorio et al. (2004) but the uncertainty is 
larger because our value includes the uncertainty of the HIPPARCOS distance measurement 
and theirs does not. 



4.2. Velocities of the GL 569 B Spectroscopic Binary 

Non-AO NIRSPEC spectra of GL 569 B do not angularly resolve Ba and Bb. A TOD- 
GOR analysis of blended spectra of GL 569 B using our suite of templates (Table 2) showed 
that the template spectrum of GL 644 G (M7) provided high correlation. The analysis also 
showed that precise measurement of the component velocities would be difficult because 
the radial velocity difference between Ba and Bb is only a few km s~^, the spectra of the 
components were similar and complex, and the spectral lines of Ba appeared to be signifi- 
cantly broadened. We expected that we would obtain the highest velocity precision by using 
angularly resolved spectra of Ba and Bb themselves as templates. To this end, we used 
NIRSPEG with AO to measure their angularly resolved spectra on UT 20 April 2003 and 
UT 24 May 2004. The small angular separation of Ba and Bb required care in the extraction 
of their spectra to minimize contamination. Figure 6 illustrates our procedure. The lower 
panel shows a sample profile of the two-dimensional spectrum in the cross-dispersion direc- 
tion from UT 24 May 2004, obtained by median filtering a REDSPEC rectified image along 
the dispersion direction; the AO plate scale is 0'.'0182 pixel"^. There is significant blending 
in the region between Ba and Bb. We modelled the profile as the sum of two Gaussians of 
equal width with separation equal to that of Ba and Bb. The upper panel of Figure 6 shows 
the extraction widths, denoted by the vertical lines and hashed regions, chosen to avoid the 
central region where the profiles are blended. The FWHM of the Gaussians in Figure 6 is 
about twice the diffraction limit at 1.6/im. The quality of the seeing and the AO correction 
varied during the May 2004 observations; the PSF in the profile shown was especially broad. 
We varied the extraction widths to suit each observation. Figure 6 shows that even in this 
case the selected regions suffer little contamination. We believe therefore that our procedure 
produced individual target spectra and hence templates useful for analysis of the blended 
spectra. 

To determine the radial velocity of Ba in orders 48 and 49 on UT 20 April 2003 and 
UT 24 may 2004, we cross-correlated spectra of Ba with our GL 644 G template, and 
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found the highest correlation for Vrotsini — 25 km s~^, if we interpret the hne broadening 
as attributable to rotation. Similarly, for Bb, using the GL 644 C template, the highest 
correlation was obtained with Vrotsini = 10 km s^^. The measured radial velocities of Ba 
and Bb in the two orders on the two dates differed by the measurement error. To use the 
individual Ba and Bb spectra as templates to analyze the blended, non-AO B component 
spectra, we averaged the measured velocities in orders 48 and 49; these values appear in 
Table 7. 

Figure 7 shows, in the heliocentric reference frame, the April, 2003 and May, 2004 AO 
resolved spectra of GL 569 Ba and Bb in orders 48 and 49 and, for comparison, a non-AO 
spectrum of GL 569 B measured on UT 22 February 2005. The spectral range displayed for 
order 48 is less than that for order 49 because we did not analyze the ~0.007 iim wide section 
shortward of 1.585 //m that suffers from terrestrial CO2 absorption (see Bender et al. 2005, 
Figure 1). It is evident that the lines in the spectrum of Ba are more broadened than those of 
Bb. Except for the effects of broadening, the spectral features in Ba and Bb appear identical 
suggesting that the spectral types must be very similar. These results are consistent with 
the finding of Zapatero Osorio et al. (2004) that, in their high-resolution, angularly resolved 
J— band spectra, also obtained using NIRSPEC, Ba and Bb appeared as spectral type M8.5 
and M9, with significant broadening characterized by v^otsini ~ 37 km s'^ and 30 km s~^, 
respectively. 

Wc analyzed the non-AO spectra with TODCOR twice using the two different AO- 
rcsolvcd templates for Ba and Bb from the two different NIRSPEC with AO runs. Wc list 
the radial velocities in Table 7. The OH emission lines in order 48 are weaker than those in 
order 49; on UT 22 December 2002 and UT 14 June 2003, the integration times were not 
long enough to expose the OH lines for a reliable dispersion solution and we do not list order 
48 velocities for those dates. We considered the difference in velocities derived on a given 
date using the two sets of AO-resolved templates as an empirical measure of the error. Thus, 
for each order, we considered the differences 

l^Ba,i(4/03) - VBa,i{^m 

in which the Vea.j's are the radial velocities of Ba measured on the ith occasion using the AO 
templates obtained in April, 2003 and May, 2004 (Table 7). The same procedure was applied 
to the velocities of Bb. The distribution of the 28 differences is centered at 0.1 ± 1.7 km s^^, 
indicating that there is no systematic velocity offset in either set of templates. The standard 
deviations for the l^Ba,i and ^Bb,i are 1.8 and 1.4 km s~^, respectively; these values are 
hsted in Table 7 as the velocity uncertainties of the Ba and Bb components. We attribute 
the larger uncertainty associated with the measurements for Ba to the greater velocity width 
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of its spectral lines. 

4.3. Dynamical Masses of the Components 

Velocity measurements of a double-lined spectroscopic binary (SB2) yield P,e,To,u!, 
asini, in physical units, and the velocity semi-amplitudes of the primary and secondary, 
Ksa and Ksb- These are related to the mass ratio, q — Msb/MBa, a, and UBa and asb, the 
semi-major axes of the primary and secondary by g = Ksa/ Ksb and asini — {aBa + o,Bb)sini. 
ttBaSini is given by 

asasini = 0.01375irB„P(l - e^)^/^ 

where asa is in gigameters, in km s~^, and P in days. The same result holds for UBb in 
terms of Ksb (Heintz 1978). 

Here, P,e,i,a{"), and To arc already known to high precision from the astrometric 
observations of GL 569 B. In combination with the HIPPARCOS distance of GL 569 A, 
a is known in physical units (Tabic 6). Our analysis of GL 569 B as an SB2 began by 
calculating the difference, in the sense, of possible velocity solutions to the data, deriving 
^Ba-i ^Bb and 7, the center-of-mass velocity, for the minimum value of ■ We did this twice, 
the first time applying the templates for Ba and Bb measured in April, 2003, and the second 
using the May, 2004 templates. The minimum reduced x? is ~1.1 for both solutions. The 
results for the velocity semi-amplitude of Bb and the center-of-mass velocity are in excellent 
agreement with each other, but the results for the velocity semi- amplitude of Ba are not. 
Using the April, 2003 templates, 

KBa = 2.90 ± 0.50, KBb = 5.34 ± 0.50, 7 = -8.47 ± 0.30 km s"^ 
and using the May, 2004 templates, 

KBa = 4.30 ± 0.70, KBb = 5.30 ± 0.50, 7 = -8.50 ± 0.30 km s~^ 

We attribute the large difference of the two KBa values to the large velocity width of 
Ba's spectrum, as also manifested (§4.2) in the scatter of the ^Ba,i-, larger than that of the 
Afib^j. Because KBb — 5.3 ± 0.5 km s'^ is well-determined, we analyzed GL 569 B as single- 
fined system on the basis of the Bb data. Figure 8 shows the velocity vs phase data for the 
April, 2003 and May, 2004 templates. The velocities of the secondary are fitted with KBb 
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and 7 as given above. Using the April, 2003 data, the mass function f{M) of the secondary 
is. 



{Mna + MBbf 



= 1.036 X 10-^is:|bP(l - e2)i Mq 

with K in km s~^ and P in days (adapted from Heintz 1978). Using Ksb = 5.3 ±0.5 km 
and values for P, e, i and Msa + Msb from Table 6, we derive that M^a = 0.105 ± 0.011 Mq. 
By subtraction from the total mass, Msb — 0.020 ± 0.013 M : . 

The value of M^a cannot be taken at face value given the similarity of the Ba and Bb 
spectra (Figure 7) and their flux ratio ~ 0.6 (Table 5). A promising explanation is that 
Ba itself is a close, angularly unresolved binary, and that GL 569 B consists of three brown 
dwarfs, two at Ba, one at Bb as previously suggested by Martin et al. (2000) and Kenworthy 
et al. (2001). The masses of the Ba brown dwarfs are probably ~0. 045 — 0.050 Mq each, and 
that of Bb, ~0.020 — 0.040 Mq. This explanation would account for the observed facts: that 
Ba is about twice as bright as Bb at H, that the total mass is ~0.12 M0, and that the spectra 
of Ba and Bb appear identical except for the greater line width of Ba. The greater line width 
of Ba could be attributable either to rotational broadening, which could disguise the presence 
of an additional component, or some combination of line broadening and Ba orbital motion. 
If the candidate binary possesses an orbital inclination distinct from and larger than that 
of the Ba— Bb system, the confusion between broadened lines and velocity shifts would be 
accentuated, particularly for objects of the similar spectral type. More AO-resolved, high 
resolution spectra of Ba are required to distinguish among these possibilities. 

Figure 9 shows GL 569 Bb on a mass-luminosity (at H) diagram; if its mass at the la- 
upper bound is 0.033 M©, the isochrones indicate a ~100 MY age, consistent with that of 
GL 569 A. We do not know the individual if— magnitudes or masses of the brown dwarfs 
proposed to comprise GL 569Ba. If their masses are ~0.045 M0 with M^^ a few tenths 
brighter than 11 mag, the isochrones indicate an age > 100 MY. The uncertainties are 
suffficiently large that it is unclear whether these mass estimates are in conflict. 

The derived center-of-mass velocity of Ba— Bb is well determined at 7 = —8.50 ± 
0.3 km s~^. It differs, however, by 3 km s~^ from the value reported by Zapatero Oso- 
rio et al. (2004), —11.5 ± 0.5 km s^^. Unlike our measurement, which is referenced to the 
stellar radial velocities of the spectral type templates, Zapatero Osorio et al. (2004) measured 
the Doppler shift of the K I doublet in their J-band spectra using the observed and labo- 
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ratory wavelengths. We believe that our value is representative of the actual value because 
it is based on all the lines in the spectra. Our measurement thus avoids the complications 
of measuring a Doppler shift of the pressure broadened and nearly saturated lines of the 
K I doublet. Moreover, for M spectral types, contamination by Mg 1 at 1.24369 fim and 
particularly by Cr I at 1.25253 /im for later types can shift the line centers of the K 1 doublet 
at 1.24356 /im and 1.25255 /xm. The fact that our measured center-of-mass velocity differs 
from the radial velocity of GL 569 A by only 0.8 ±0.9 km s~^ also suggests that it is close to 
the actual value. The ~3 km s"^ difference between our value and that of Zapatero Osorio et 
al. (2004) affects only their absolute radial velocities. The differences in the radial velocities 
of Ba and Bb measured by Zapatero Osorio et al. (2004) (their Table 4 and Figure 8) agree 
well with the differences derived from our values in Table 7, as shown in our Figure 10. 



5. Discussion 

5.1. The GL 569 Quadruple? 

Our work suggests that GL 569 may be a quadruple consisting of an M2.5 V star 
(GL 569 A) accompanied by three nearly equal mass brown dwarfs and distributed as follows: 

1) The wide binary, GL 569 A — GL 569 B, at separation ~5" (~50 AU) with total mass 
(0.35 ± 0.03 + 0.125 ± 0.007) = 0.48 ± 0.03 Mq. If it is oriented face-on, its orbital period is 
~500 yrs. 

2) The binary GL 569 Ba - GL 569 Bb which has semi-major axis 90 mas (0.89 AU), 
period 864 days, and total mass 0.125 ± 0.007 Mq. Bb is a single brown dwarf with mass 
0.020 ±0.013 Mq. 

3) Ba is a candidate unresolved inner brown dwarf binary with total mass 0.105 ± 0.011 Mq. 

Designating the two components of Ba, the inner binary, as Baa and Bab, we consider 
the limits that can be placed on their orbit empirically by our spectroscopic and imaging 
observations, and theoretically by the requirement that Ba and Bb form a stable triple. By 
simulating the inner binary as composed of two equal brightness brown dwarfs of late M 
spectral type, rotationally broadened to 10 km s~^, we found that we could measure their 
individual velocities reliably to ~5 km s~^. We measured no significant velocity differences 
in the two AO spectra of Ba. We conclude that the observed velocity difference of the two 
components was < 5 km s~^ at the times of the observations if the components had ~ 10 
km s^^ rotational broadening. With Mfot = 0.105 Mq, this implies that the semi- major axis 
of the inner binary satisfies ainner ^ 3.8 sin^ij„„er AU, where imner is the incUnation of the 
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inner binary. 

The FWHM of the H-band PSF was 41 mas in the images obtained on UT 25 February 
2005. We analyzed images of model binaries with equal brightness components convolved 
with the PSF and the same signal to noise as the H-band image in Figure 4. We were able 
to identify the two components at separations as small as ~ 10 mas. The actual situation is 

complicated by the presence of the third component, Bb, at apparent separation ~ 85 mas 
(Table 5), 2 x FWHM of the PSF. It is apparent in Figure 4 that the broad wings of the 
PSF extend from Bb to Ba. We estimate that Baa and Bab would have to lie at apparent 
separation > 20 mas for their presence to be recognizable in our images. This limit, combined 
with that from the spectroscopic observations, restricts iinner to < ±13.6°. 

Theoretical study of the stability of triples has an extensive literature and continues to 
be an interesting area of research as new environments are identified (e.g., the formation 
of planets in stellar binaries, Holman & Weigert 1999) and the meaning of "stability" is 
sharpened (e.g.. Ford et al 2000). Eggleton & Kiseleva (1995) presented an empirical stability 
criterion in terms of y"**", the ratio of the periastron distance of the outer binary to the 
apastron distance of the inner binary. 



yrmin QoMterl,! ^outer) 

^inner{^ ~l~ ^inner) 

When values of y™*" larger than a critical value yj"*" are realized, the triple is stable. The 
criterion is applicable to triples with various inclinations, eccentricities, and relative phases 
of the inner and outer binaries. For a triple with equal mass components (see Eggleton & 
Kiseleva 1995, eqn. 2), Y^'"" ~ 5. With a^ter = 0.90 AU and e^ter = 0.31 (Table 6), 

0'inner{\ ~l~ dinner) — 0.13 AU; dinner — 0.13 AU if G-inner — 0, and 0.07 AU if G-inner — 0.9. In 

the latter case, the apastron of the inner binary is 0.13 AU. We conclude that the largest 
value of dinner for Stability, 0.13 AU, is consistent with our observational limits if iinner lies 
in the range ±11°. 



5.2. The Lithium Test 

The Li spectral line at 6708 A is a sensitive mass and age indicator because brown 
dwarfs with mass > 0.06 destroy their Li within a few hundred Myr (Magazzu et al. 
1993; Chabrier & Baraffe 2000). The masses we have estimated for Baa, Bab, and Bb are 
at this boundary. In spectroscopic observations that did not resolve the GL 569 B system, 
Magazzii et al. (1993) did not detect the Li line, suggesting that the brown dwarf masses 
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are greater than the values we infer. Detection of the Li feature in GL 569 B may have 
been frustrated by the large line widths revealed by our analysis of the IR spectra and the 
possible triple nature of the system. 

6. Summary and Suggestions for Future Work 

There are two main results of our work: 

1) The previous suggestions that Ba could be a binary brown dwarf were based on the 
fact that Ba is ~ 1.6 times brighter than Bb in the near IR and that the near IR colors of 
the two are similar (Martin et al. 2000; Kenworthy et al. 2001). Our work confirms these 
observations and strengthens the case for the binarity of B a in three respects: i) the finding 
that Mfia = 0.105 ±0.011 Mq and Ba cannot be a star of this mass, ii) the similar spectra of 
Ba and Bb indicate that Ba and Bb contain similar brown dwarfs, and Hi) the greater width 
of the lines of Ba, ~ 25 km s^^, than that of Bb, ~ 10 km s~^, has a natural explanation if 
Ba harbors two brown dwarfs. 

These arguments and the total mass of the system Mtot = 0.125 ±0.007 M© also suggest 
that the masses of the three components Baa, Bab, and Bb, are approximately equal and in 
the range 0.04 to 0.05 Mq. 

2) The location of GL 569 A with respect to theoretical isochrones and in the Pleiades 
CMD, and the location of Bb with respect to theoretical isochrones (and also of Baa and Bab 
if they have approximately equal mass) indicate that the age the GL 569 system is 100-125 
Myr. 

The GL 569 quadruple is well suited to improve our understanding of objects at the low 
mass end of the IMF. We suggest the following observations to confirm and advance beyond 
the results we have presented. 

A) Detection of the Li line in the spectra of GL 569 Ba and Bb both would confirm the 
mass of GL 569 Bb and that Ba must be composed of two brown dwarfs. The spectroscopic 
measurement should resolve Ba and Bb to avoid the complications of interpreting the complex 
spectrum of three unresolved brown dwarfs. This would require AO-assisted spectroscopic 
observations in the visible, a capability not yet available. 

B) Additional high angular resolution IR spectroscopy of Bb will improve the precision of 
the velocity semi- amplitude Ksb, and hence of Mbq and Mb&. 

C) High angular resolution IR spectroscopy of Ba may distinguish the two brown dwarfs Baa 
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and Bab, and hence determine their parameters as a double-hned spectroscopic binary. 

D) At separation ~7 mas, the inner binary is resolvable by the present generation of IR 
interferometers and is not far below their present sensitivity limits. The VB parameters of 
the inner binary, combined with the SB2 parameters would determine the individual masses 
Msaa and Msab ^nd reveal whether the orbits in the hierarchical triple are coplanar. 

Of the observations suggested only the AO-resolved IR spectroscopic observations B and C 
are possible now. We plan to continue with these and look forward to the technical advances 
that will enable the interferometric and visible light spectroscopic studies. 
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Table 1. Log of Observations 



UT Date 


MJD 


Instrument 


Mode 


Targets 


2002 Jul 18 


52473.3 


NIRSPEC 


non- 


■AO 


GL569B 


2002 Dec 22 


52630.7 


NIRSPEC 


non- 


■AO 


GL569B 


2003 Mar 25 


52723.5 


NIRSPEC 


non- 


■AO 


GL569B 


2003 Apr 20 


52749.4 


NIRSPEC 


AO 




GL569A,GL569Bab 


2003 Jun 14 


52804.2 


NIRSPEC 


non- 


■AO 


GL569A,GL569B 


2004 Jan 26 


53030.6 


NIRSPEC 


non- 


■AO 


GL569B 


2004 Apr 29 


53124.4 


NIRSPEC 


non- 


■AO 


GL569B 


2004 May 24 


53149.3 


NIRSPEC 


AO 




GL569Bab 


2004 May 25 


53150.3 


NIRSPEC 


AO 




GL569A 


2004 Dec 24 


53363.7 


NIRC2 


AO 




GL569Bab 


2004 Dec 26 


53365.7 


NIRSPEC 


non- 


■AO 


GL569B 


2005 Feb 22 


53423.5 


NIRSPEC 


non- 


■AO 


GL569B 


2005 Fel) 25 


53126.5 


NIRC2 


AO 




GL569A.GL569Ba,b 
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Table 2. Revised Template Library Radial Velocities 



Star 


Sp. Type 


Obs 


. Date 


Vrad (km s ^) 


(jLioA 


Mi. 5 


zOOO 


Jun 10 


11 1 

11. el 


GL436 


M2.5 


2001 


Jan 5 


9.61 


GL752A 


M3 


2000 


Jun 9 


35.88 


GL213 


M4 


2000 


Jan 8 


105.96 


GL402 


M4 


2001 


Jan 5 


-1.04 


GL669B 


M4.5 


2000 


Jun 10 


-34.90 


GL406 


M5.5 


2001 


Jan 5 


18.61 


LHS292'^ 


M6.5 


2001 


Jan 5 


1.35 


GL644C 


M7 


2001 


Feb 2 


15.21 


LHS2351 


M7 


2001 


Feb 2 


0.72 


LHS2065'' 


M9 


2001 


Feb 2 


7.02 


2MASS 0208+25 


LI 


2002 


Dec 15 


20.79 



^Reported to be a spectroscopic binary; ''X-ray flare star 
(Guenther & Wuchterl 2003) 





Table 3. 


Astrometry of GL 569 AB 


Date 


R.A. X cos(Dec) (") 


Dec (") Reference 


1985.6 


1.52 ±0.15 


4.84 ± 0.15 Forrest, Skrutskie, & Shure (1988) 


1999.7 


2.11 ±0.08 


4.53 ± 0.04 Martin et al. (2000) 


2001.1 


2.45 ±0.22 


4.23 ± 0.13 Lane et al. (2001) 


2005.2 


2.72 ±0.14 


1.18 ±0.21 this work 
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Table 4. Radial Velocity Measurements of GL 569 A 



MJD 


Vrad (km s ^) 


52749.4 


-7.1±1.0 


52804.3 


-8.1±1.0 


53150.3 


-7.8±1.0 


Average 


-7.7±0.6 



Table 5. Photometry and Astrometry of GL 569 B 



M.li) 


Ex]). (s) 




llDl. 




Kbi, 


Amn 


AuiK 


Sep. [") 


1>. A. (-) 


53363.7 


100 










0.61 ±0.03 




0.0885 ± 0.0048 


111.3 ± 1.2 


53426.5 


22 


10.43 ± 0.04 


11.04 it 0.05 


9.86 ± 0.05 


10.39 ± 0.06 


0.57 ±0.04 


0.61 ± 0.04 


0.0798 ± 0.0029 


133.7 ±0.7 



Table 6. Orbital Parameters of GL 569 Bab 



Parameter Visual Binary Orbit 

Zapatero Osorio et al. (2004) This Work 



P (days) 


876 ±9 


863.7±4.2 


T (MJD, days) 


51, 822 ±3 


51820.9±2.6 


e 


0.32 ±0.01 


0.312±0.007 


a (mas) 




90.4±0.7 


a (AU) 


0.90 ±0.0P 


0.89 ±0.02*^ 


i (dcg) 


34 ±2 


32.4±1.3 


u (deg) 


257 ±2 


256.7±1.7 


Q (deg) 


321.5 ±2.0 


321.3±2.2 


Mtotal (Mo) 


0.125 ±0.005^^ 


0.125 ±0.007^ 



'^Docs not include 1.6% parallax uncertainty; '^does include 
1.6% parallax uncertainty 
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Table 7. Radial Velocity Measurements of GL 569 Bab 



MJD phase V^ad (km s ^) 

GL 569 Ba GL 569 Bb 

Order 49 Order 48 Order 49 Order 48 



GL644C Templates 



52749.4 0.072 -7.92±0.90 -9.98±0.90 

53149.3 0.535 -8.50±0.90 -8.38±0.90 



April 20, 2003 Templates 



52473.3 


0.753 


-10.89±1.80 


-11.60±1.80 


-4.71±1.40 


-3.19±1.40 


52630.7 


0.935 


-11.60±1.80 




-3.65±1.40 




52723.5 


0.042 


-9.46±1.80 


-7.12±1.80 


-8.07±1.40 


-7.74±1.40 


52804.2 


0.136 


-7.21±1.80 




-14.15±1.40 




53030.6 


0.397 


-7.52±1.80 


-4.92±1.80 


-12.72±1.40 


-11.74±1.40 


53124.4 


0.506 


-9.07±1.80 


-8.54±1.80 


-9.56±1.40 


-7.29±1.40 


53365.7 


0.785 


-10.62±1.80 


-12.31±1.80 


-4.98±1.40 


-3.09±1.40 


53423.5 


0.852 


-10.75±1.80 


-11.71±1.80 


-3.99±1.40 


-1.98±1.40 


May 24, 2004 Templates 


52473.3 


0.753 


-12.30±1.80 


-10.93±1.80 


-4.57±1.40 


-3.38±1.40 


52630.7 


0.935 


-12.12±1.80 




-3.96±1.40 




52723.5 


0.042 


-6.92±1.80 


-4.54±1.80 


-10.28±1.40 


-9.98±1.40 


52804.2 


0.136 


-6.29±1.80 




-13.21±1.40 




53030.6 


0.397 


-5.42±1.80 


-4.49±1.80 


-13.03±1.40 


-11.28±1.40 


53124.4 


0.506 


-6.96±1.80 


-5.34±1.80 


-11.64±1.40 


-11.14±1.40 


53365.7 


0.785 


-12.90±1.80 


-12.66±1.80 


-4.40±1.40 


-3.75±1.40 


53423.5 


0.852 


-12.76±1.80 


-11.82±1.80 


-3.41±1.40 


-2.70±1.40 



-20- 



6 



o 

Q) 
(fi 
O 

o 



Q 



5 
4 
3 
2 
1 



GL569 B 

July 1985 -f 

♦ Aug 1999 
Feb 2005 Feb 2001 





-1 



4 3 2 
R.A. (arcsec) 



GL569 A 
1 -1 



Fig. 1. — The position of GL 569 B with respect to A using our measurement and the earher 
measurements hsted in Table 3. A and B form a common proper motion system (see §3); 
the figure shows orbital motion of B with respect to A. North is up and east is to the left. 
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Fig. 2. — Absolute V and H magnitudes and mass of GL 569 A compared to theoretical 
isochrones calculated by Baraffe et al. (1998). 
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Fig. 3. — Absolute K mag and V— K color of GL 569 A (filled circle) on a color-magnitude 
diagram of the Pleiades. The Pleiades data is the same as in Luhman et al. (2005) Figure 
1, and is for an adopted extinction Ay = 0.12 mag and distance 133 pc for the Pleiades. 
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Fig. 4. — NIRC-2 adaptive optics images of GL 569 Ba and Bb at H and K obtained on 
UT 25 February 2005. North is up and east is to the left. The images show a subsection 
of the frame; GL569 A was imaged in the same field, but is not shown. The first contour 
represents 1000 counts, ~ 5(T over the background noise, and each subsequent contour is an 
additional 1000 counts. Ba is located northwest of Bb and is clearly the brighter component; 
the Bb/Ba flux ratios are 0.57 ± 0.04 and 0.61 ± 0.04 in H and K respectively. 
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Fig. 5. — The astromctric orbit of GL 569 Bb with respect to Ba at the coordinate origin, 
using our measurements from UT 24 December 2004 and UT 25 February 2005 together 
with the earher measurements given in Zapatero Osorio et aL (2004). The best fit orbit 
corresponds to our updated parameters in Table 6. 
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Fig. 6. — A sample profile in the cross-dispersion direction of the AO angularly resolved 
Order 49 GL569 Ba and Bb two-dimensional spectrum from UT 24 May 2004. The open 
circles (bottom panel) show the measured spatial profile of the resolved spectra; the solid line 
is a least squares fit of two Gaussians with equal widths (see text) . The dotted and dashed 
curves (top panel) are the individual Gaussian profiles. The vertical lines and hashed areas 
indicate the spatial regions used to extract each spectrum. 
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Fig. 7. — Angularly resolved spectra of GL 569 Ba and Bb measured on UT 20 April 2003 
and UT 24 May 2004, compared with a spectrum measured on UT 22 February 2005 in which 
Ba and Bb are unresolved. NIRSPEC spectral orders 48 (left) and 49 (right) are shown. The 
wavelength span of the order 48 spectra is about 2/3 that of the order 49 spectra because the 
~0.007//m wide section shortward of 1.585 //m is contaminated by a terrestrial CO2 band 
and not shown. The spectra are normalized to 1 in the continuum and are plotted on a to 
1 flux scale with the Ba spectra in the top two panels offset by +1 in flux. 
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Fig. 8. — Velocity vs orbital phase measurements for GL 569 Ba (open triangles) and Bb 
(filled triangles) using the UT 20 April 2003 templates (left) and UT 24 May 2004 templates 
(right). Order 48 and 49 data are plotted as downward and upward pointing triangles, 
respectively. The measurements obtained with AO, averaged over orders 48 and 49, are 
indicated with open (Ba) and filled (Bb) circles. The dashed fine is the fitted center-of-mass 
velocity, —8.5 km s~^, and the solid line is the model fit with Kbi — 5.3 km s~^ and P, e, i, 
and cu given in Table 6. 
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Fig. 9.— GL 569 Bb with M^^ = 11.08 ± 0.05 mag and mass 0.020 ± 0.013 M© compared 
with theoretical isochrones calculated by Baraffe et al. (2003). Its location is consistent with 
an age of ~ 100 Myr. 
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Fig. 10. — Differences of the GL 569 Ba and Bb velocities vs orbital phase using our measure- 
ments in Table 7 with the UT 20 April 2003 templates for: order 48 (downward triangles), 
order 49 (upward triangles), and AO (circles). Our data are consistent with the values 
(squares) reported by Zapatero Osorio et al. (2004). 



